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Abstract 

The q- Dobinski formula may be interpreted as the average of pow- 
ers of random variable X q with the q-Poisson distribution. 

Forty years ago Rota G. C. 1 proved the exponential generating function 
for Bell numbers B n to be of the form 

oo n 

J2^(B n )=exp(e*-l) (1) 

n=0 

using the linear functional L such that 

L(X^) = 1, n > (2) 

Then Bell numbers (see: formula (4) in pQ) are defined by 

L(X n ) =B n , n > (3) 

The above formula is exactly the Dobinski formula [2| if L is interpreted as 
the average functional for the random variable X with the Poisson distri- 
bution with L(X) = 1. It is Blissard calculus inspired umbral formula 
. Recently an interest to Stirling numbers and consequently to Bell num- 
bers was revived among " g-coherent states physicists" EJ E] • Namely the 
expectation value with respect to coherent state I7 > with [7] = 1 of the 



1 



n-th power of the number of quanta operator is "just" the n-th Bell number 
B n and the explicit formula for this expectation number of quanta is "just" 
Dobinski formula ,3J. The same is with the g-coherent states case jH] i.e. the 
expectation value with respect to g-coherent state I7 > with I7I = 1 of the 
n-th power of the number operator is the re-th g-Bell number [SJ |3] defined 

as the sum of g-Stirling numbers due to Carlitz as in Note 

there then that for standard Gauss q-extension x q of number x we have 

n 

<=E{I}„4 w 

Hence the expectation value with respect to g-coherent state I7 > with 
I7I = 1 of the n-th power of the number operator is exactly the popular q- 
Dobinski formula . It can be given via (3) Blissard calculus inspired umbral 
formula form and may be treated as definition of B n (q) 

L q (X%) = B n (q), n > 0. (5) 

due to the fact that linear functional L q interpreted as the average functional 
for the random variable X q with the g-Poisson distribution with L q {X q ) = 1 
satisfies 

L q (Xf) = 1, n > 0. (6) 

We arrive to this simple conclusion using Jackson derivative difference op- 
erator in place of D = d/dx in q =1 case and the power series generating 
function G(t) for g-Poisson probability distribution: 

\ k 

Pk = [exp.X]- 1 — G(t) = Y,Pkt\ (7) 

K q . 



n>0 



Pn = [-*— T^]t=o, [d q G(t)] t=1 = lforX = 1. (8) 



n q . 



There are many g-extensions of Stirling numbers according to their weighted 
counting interpretation. For example w(tt) = q^^^^^w^) = q mv ( n ) from 
[7] gives after being summed over the set of /c-block partitions the Carlitz 
^-Stirling numbers or w(ir) = q nin ( 7r ) from jS] gives rise to Carlitz-Gould 
q- Stirling numbers after being summed over the set of A:-block partitions 
or with w(tt) = q 1 ^ in jH] - we arrive at another combinatorial interpreta- 
tion of ^-extended Stirling numbers. g-Stirling numbers much different from 
Carlitz g-ones were introduced in the reference JHJ from where one infers 
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jllj the cigl-ansdog of (5) . Let II denotes the lattice of all partitions of the 
set {0, 1, .., n — lj.Let ir £ LT be represented by blocks ir = {B Q , B\, ...Bi, ...} 
, where B Q is the block containing zero: € B . The weight adapted by 
Cigler defines weighted partitions' counting according to the content of B Q . 

Namely w(tt) = , cigl(ir) = £ ZeBo Z, £.eA n , fc 9 cigl{n) = {l} q there- 

fore X^en q c%91 ^ = B n (q). Here A n ^ stays for subfamily of all A;-block 
partitions. With the above relations one has defined the cigl-q-Stirling and 
the cigl-q-Jiell numbers. The cigl-q-Stir\mg numbers of the second kind 
are expressed in terms of g-binomial coefficients and q = 1 Stirling num- 
bers of the second kind [EH • These are new g-Stirling numbers. The 
corresponding cigl-q-Bel\ numbers recently have been equivalently defined 
via czg/-g-Dobinski formula ^T] L(X™) = B n (q), n > 0, X™ = 
X(X + q - 1)...(X - 1 + q n ~ l ) interpreted as the average of this specific 
n — th cigl-q-power random variable X™ with the q = 1 Poisson distribution 
such that L(X) = 1. To this end note that in , JH] a family of the so 
called ip- Poisson processes was introduced. The corresponding choice of the 
function sequence ip leads to the g-Poisson process. Accordingly the exten- 
sion of Dobinski formula with its elementary essential content and context 
to general case of ip- umbral instead g-umbral calculi case only - is automatic 
in view of an experience from ^2] , [13] (see corresponding earlier references 
there and necessary definitions). At first what you do is to replace index q 
by ip in formulas (3), (4),..., (8) . Then you have got started problems with 
not easy combinatorial interpretation if at all and... etc. ^-Stirling numbers 
and -0-Bell numbers are being then defined by (4) and (3) correspondingly 
with q replaced by ip. We get used to write these extensions in mnemonic 
convenient upside down notation ^2] , 

ip n = n^,x^ = ip(x) = ip x ,n^ = n^(n - 1)^,1, n > 0, (9) 

x-p = x^(x - l)$(x - 2)^...(x - k + 1)^ (10) 

x^(x - l)^...(x —k + l)^ = ip{x)ip{x - l)...ip(x - k - 1). (11) 

You may consult for further development and use of this notation |12j , |13j 
and references therein. 
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